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By choosing two typical input parameter points in the minimal supergravity 
\ (mSUGRA) model and using the QCD factorization (QCDF) approach, we studied 

CN ' the supersymmetric effects to the CP violation of the two-body charmless hadronic 

■ B meson decays. We found that though the SUSY contributions can give large 

■ corrections to the CP asymmetries for some decay channels, they could not be dis- 
tinguished experimentally from the SM values because of the large theoretical errors 

£S) \ dominated by calculating the annihilation contributions in the QCDF approach. 
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As one of the main goals of B experiments, measurements of CP asymmetries in various 
t— i \ B decay processes are at the center of attention. Since in the standard model (SM) all 
the measured CP asymmetries have to be consistently explained by the complex phase 
5 in the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix, precise measurements of 
CP asymmetries at the current B factories can be used to test the CKM mechanism and 
C^; sequentially the SM[i|,0 

Experimentally, Belle collaboration has given the first experimental evidence for the 
direct CP violation in the decay mode B — > 7r + 7r~[3]. And presently large direct CP 
violation in neutral B meson decay B° — > K + n~ has been measured by both BaBar[4j 
and Belle jEJ collaborations, but no CP violation has been found for the charged B meson 
X ■ decay B+ -> K + tt°. The world averages given by the Heavy Flavor Averaging Group 
£ i (HFAG) are 

A CP (B° -> K + ir~) = -0.115 ±0.018, 

A CP (B + -> K + ir°) = -0.04 ±0.04. (1) 

In the SM and with the QCD factorization (QCDF) approach, the CP-violating asym- 
metry for these two channels should be naturally very similar in size and positive in sign 
for a set of favored hadronic parameters @. In the perturbative QCD (PQCD) factor- 
ization approach^, however, the SM predictions are [H Acp(B° — > K + n~) w —0.18 and 
Acp(B° — > K + n~) w —0.15: the first prediction agrees well with the measured value as 
given in Eq. (fl|). but the second one shows a large deviation with the data. 

For the measurement of the time-dependent CP asymmetries, great efforts have been 
made by BaBar, Belle and other collaborations. The recently reported measurements of 
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time-dependent CP asymmetries in B — > 0X S decays and other s-penguin modes lead to 



sin{2l3) ave = 0.50 ± 0.06, (2) 

which is a 2.6<r deviation from that for B — > J/ipK s and other charmonium modes, 

sin(2(3) ave = 0.687 ± 0.032 . (3) 

Therefore, something must be there to account for the large deviations. Since the B — > 
<f)K s and other s-penguin modes are induced at the one loop level, while B — > J/ ipK s and 
other charmonium modes are governed by tree level decays, it is tempting to expect that 
the new physics (NP) contributions are more significant to the former and do not show 
up clearly to the latter. 

Recently, many works about various NP models have been done in the B meson system 
In our previous papers we have studied the minimal supergravity (mSUGRA) 
model by calculating the branching ratios of B — > M\Mi (Mi stands for the light pseudo- 
scalar (P) or vector (V) mesons) decays, and have found some interesting results. Based on 
the previous works |l0J], here we will try to calculate how much can the SUSY contributions 
in the mSUGRA model affect the CP asymmetries of the B — > M\M 2 decays by using 
the same input parameters as listed in the Appendix of Ref . ^(J . 

In the mSUGRA model ill], only four continuous free parameters and an unknown 
sign are left, say tan/?, mi, m , A Q and sign(n). The new effects on the rare B meson 
decays from this model will manifest themselves through the corrections to the Wilson 
coefficients of the same operators involved in the SM calculation. In Ref. ^(j, we have 
found that the SUSY contributions to the Wilson coefficients Ck{k = 3 ~ 6) are always 
small and can be neglected safely. But for the Wilson coefficients C 7l {m b ) and C$ g (mb), 
the SUSY contributions can be rather large and even make them have an opposite sign 
with the SM ones. 

Similar to Ref. considering the constraints coming from the well measured 

B — ► A s 7 decays and so on, we take two typical parameter points, say Case A and 
B as listed in the following tabular. For Case C in Ref. [10], the Wilson coefficient 
C 7 ^(m b ) in the mSUGRA model is also SM-like (negative) as that in Case A. Hence in 
our calculations we will only consider Case A and B. The numerical values of the ratio 
Rj = C 7l {m b ) I (m&) for the two cases have also been given in the tabular. It is easy 
to see that the Wilson coefficient CV 7 (m&) in the mSUGRA model is SM-like (negative) 
for Case-A, but nonstandard (positive) for Case-B. 



Case 


m mi 


A, 


tan/3 Sign[n] 


R 7 


A 


300 300 





2 


1.10 


B 


369 150 


-400 


40 + 


-0.93 



Just take a look at the two typical sets of SUSY parameters, both of them are chosen 
to be real for the constraint condition coming from the measured electric dipole moment 
(EDM) of neutron and electron. Therefore no new weak phase has been introduced and 
the mechanism of CP violation in the mSUGRA model is the same as that in the SM. 
But when we calculate the amplitudes of B decays by using the QCDF approach 0, 
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the real and imaginary part of the factorized coefficient ai(MiM 2 ) or aj(M 1 M 2 ) (See 
Ref. [H for detailed definitions) will be modified by the SUSY contributions and hence 
the CP-violating parameter will be affected in this model. 

We begin with a discuss of the CP asymmetries in B — > PP decays. Firstly, for those 
channels which have only direct CP violations, one can define 

T(B -> /) + T(B -> /) <■ ' 



As we have discussed in Ref. |10j , the potential SUSY contributions are mainly embod- 
ied in a\(MiM2) and a^ ew (MiM 2 ). Therefore, we naturally expect a large new physics 
corrections to those penguin dominated B meson decays. According to our calculations, 
things are absolutely so. For the tree-dominated decay channels, say B ± — > 7r ± (7r°, 77")), 
the SUSY corrections are less than 10% for both Case A and B. However for other channels 
which are penguin-dominated, we can see from Table UJ (a) in Case A, all the corrections 
are still small and less than 10% but in Case B the corrections are large and the largest 
one is about —35% for B — ► -k K^. (b)Though the corrections in Case B are large for 
those channels, experimentally they can not be distinguished from the SM values for they 
are all of the same order of magnitude as SM values, (c) For the very interesting decay 
B — > 7r ± A =F , Aqp in both the SM and mSUGRA model have a different sign with the 
data in Eq. (JTJ. If we consider the large uncertainties from the annihilation contributions, 
A CP is 

r 4.5±|:|, sm 

Acp^RT) = } 4.9±?o 5 2 , Case - A (5) 
1 3.0tjj;?. Case - B 

and can only be consistent with the data within about three standard deviations. However, 
in Refs. [12|, by using the PQCD approach or considering the final state interaction (FSI) 
modifications to short-distance (SD) predictions in QCDF approach, the authors found 
the consistency between the theory and the experiments for B — > -n^K T decay. For more 
details, one can see these references. 



TABLE I: Theoretical results in percent for those penguin dominated and having only direct CP 
violations Aqp channels of B — > PP in the SM and mSUGRA model, where superscripts "f+a" 
and "f denote the predictions with or without weak annihilation contributions considered. 



Decay 


SM 


Case A 


Case B 


Modes 


a! Af+a 
■^CP ^CP 




•^CP 


B d -» 7T ± K =F 


5.36 4.54 


5.75 


3.47 


B ± 


8.28 7.45 


8.70 


6.09 




1.00 0.91 


1.03 


0.81 


B ± -» ^±7/ 


-13.7 -12.8 


-14.1 


-11.5 


B ± -» K^rj 


2.38 2.04 


2.46 


1.94 


B ± _» K±K° 


-26.8 -24.0 


-27.6 


-22.0 
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Secondly, for the CP asymmetries of the left B — > PP channels, they are time- 
dependent and can be described by 



A C p(t) 



T(B°(t) -> f) - T(B°(t) -> f) 

r(s°(f) - /) + r(s°(f) - /) 

5/ sin(AMt) - C/ cos(AMt) (6) 



where Cf and Sf represent the direct and the mixing CP asymmetry, respectively, and 
they are given by 

1 - |A CP | 2 _ 2Im(X CP ) 
L f~TT\\ — 12' b f ~ T^TTl — 12- 



The parameter Acp is defined by 



V t tV td A(B%0)^f) 
V tb V t * d A(B°(0)^fY 



Through the numerical calculations, we found that for those neutral B d — > PP decays, 
the SUSY contributions to the two CP asymmetric parameters Cf and Sf are generally 
small or moderate in both Case A and B. The largest correction is only about 26% to Cf 
of Bd — > KgiT and make it vary from 3.8% to 2.8%. Obviously, so small correction can 
not be distinguished experimentally and can be masked easily by the uncertainty coming 
from the annihilation contributions . 



TABLE II: The same as Table d but for some B — > PV channels. 



Decay 


SM 


Case A 


Case B 


Modes 




A t+a 


At 


At 


B ± - 




2.01 


1.71 


2.19 


1.27 


B ± - 


-» n°K*± 


9.26 


8.56 


9.81 


6.32 


B±- 


-» if^ 


3.19 


2.98 


3.21 


3.08 


B±- 


-» iCV 


-27.5 


-19.9 


-25.1 


-64.5 


B ± - 




-13.5 


-11.6 


-13.1 


-17.2 


B±- 




0.12 


0.16 


0.12 


0.19 


B ± - 




-6.61 


-5.90 


-6.43 


-8.21 


B ± - 


-> 


2.12 


1.75 


2.27 


1.44 


B ± - 


-> i^ic 


-52.7 


-46.0 


-56.7 


-34.4 


B ± - 




-2.55 


-2.92 


-2.43 


-4.03 


B d - 




-0.32 


-0.37 


-0.28 


-0.26 


B d - 




-15.4 


-12.9 


-18.1 


-7.34 


B d - 




-4.51 


-3.25 


-4.20 


-8.18 


B d - 




4.93 


4.45 


4.97 


4.69 


B d - 




-11.2 


-8.15 


-10.4 


-23.8 
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TABLE III: The same as Table [fl but for CP asymmetries parameter Cf and Sf of some 
Bd — > PV channels. 



Decay Modes 




SM 


Case A 


Case B 


B d - ttV 




-4.47 -3.03 
-32.1 -35.8 


-2.94 
-35.2 


-15.8 
-8.92 


B d -» TT°UJ 


C/ 
5/ 


74.4 94.6 
-62.5 32.3 


73.0 
-56.8 


61.6 
-64.9 


B d -> W° 


5/ 


6.28 1.94 
-29.3 -32.7 


-4.51 
-52.4 


47.6 
87.9 


£ d - 


5/ 


46.5 37.3 
-69.5 -63.3 


44.8 
-73.6 


58.1 
-24.6 




C/ 
5"/ 


8.36 5.48 
-41.5 -38.0 


7.14 
-39.4 


16.9 
-55.6 


B d -» 7?' (J 


C/ 
^/ 


-18.4 -17.9 
-28.9 -25.5 


-19.3 
-27.0 


-11.5 
-42.9 




Cf 
Sf 


-9.15 -7.91 
-62.1 -63.8 


-8.85 
-62.6 


-12.3 
-57.1 




Cf 
Sf 


9.71 8.26 
-89.8 -87.9 


9.36 
-89.5 


13.5 
-93.9 


B d -» 


Cf 
Sf 


-2.12 -1.85 
76.9 76.6 


-2.27 
76.9 


-1.44 
76.9 



We now take a look at B — > PV mode. Different from B — > PP mode, the SUSY 
correction to this mode can interfere with the SM counterparts constructively or destruc- 
tively [13]. Through the numerical results, we found the new physics corrections to the 
CP violations in Case B are large for most channels. For those channels having only direct 
CP violation, the penguin dominated ones are affected a lot. As shown in Table |Hj in 
Case B the largest SUSY corrections is about 135% for B — > K*^r\ channel and make 
the size of its CP violation increased from 28% to 65%. Such large corrections may be 
measured experimentally. 

For other neutral Bd — > PV decays, they have both direct CP violation Cf and mixing 
CP asymmetry Sf. For the CP-violating parameters of B — > (ir, rj^)(j) decays, since 
these channels receive no contributions from electromagnetic or chromo-magnetic penguin 
operators where the SUSY contributions are entered, they remain almost unchanged in 
the mSUGRA model. For other channels as shown in Table 11111 the direct CP violations 
of most channels in Case B are greatly affected by the SUSY corrections. The largest 
corrections even reach a factor of 7 for the decay B° — > rjp , about 253% increase for 
B° — > tt°p° and 100% enhancement for B° — > r\uo. As to the indirect CP violation Sf in 
Case B, the SUSY corrections are also large. For Bd — > r/p decay, for example, the sign 
of its indirect CP violation has been changed and the size also increased by a factor of 
3. But for the very interesting channels B — > (j)K s , the SUSY contributions make little 



6 



TABLE IV: The same as Table but for some B — > VV channels. 



A dir {B 


-» VV) 


SM 


Case A 


Case B 


A CP (B° - 




17.8 


19.3 


10.5 


Acp(B° - 


> K*°p°) 


-21.1 


-23.4 


-11.9 


A CP {B- - 


-> A'*-/j°) 


18.7 


19.6 


13.8 


A CP (B- - 


■» K*V) 


1.57 


1.63 


1.25 


C f {B° K*°K*°) 


-28.2 


-29.4 


-21.6 


A C p{B- - 


-» K*~K*°) 


28.2 


29.4 


21.6 


A CP (B° - 




12.8 


13.4 


9.68 


A C p(B- - 




32.5 


34.3 


22.4 


Acp(B° - 


► if *V) 


1.75 


1.81 


1.38 


A C p(B- - 


+ K*-0) 


1.75 


1.81 


1.38 



effects and hence the U (J)K S " puzzle as mentioned above still cannot be solved here. 

At last, let us talk about B — > VV mode in brief. For B — ► VV, we first give two 
remarks: (a)from Eq.fjHJ) one can see that the parameter Xcp for B — > VV mode is helicity- 
dependent since three decay amplitudes with A = (0, ±1) are different. It follows that 
theoretically the indirect CP violation Sf can be given only when the two vector mesons 
in the final state have a certain helicity. Hence in our calculation we only give the direct 
CP violation (Acp or Cj) of B — > VV in Table ITVl (b)the annihilation amplitude in the 
VV case does not gain a chiral enhancement of order M^/{m q m}j) as that in B —>■ PP 
and PV modes. Therefore, it is truly power suppressed in heavy quark limit and we have 
ignored such contributions in our calculation. 

Similar to B — > PP, only in the QCD penguin-dominated channels can their CP 
violation be affected a lot. From Table I1VI in Case B the mSUGRA predictions of CP 
violations for all the listed channels become smaller than corresponding SM predictions. 
As to B° -> K*°p°, the SUSY contributions in Case B can even make its CP-violatine 
parameter decreased by about 44%. However, these SUSY contributions are of the same 
order of magnitude as SM values and therefore almost impossible to be distinguished 
experimentally. 

To conclude, we have computed the CP asymmetries for two-body charmless B — > 
M 1 M 2 decays in the mSUGRA model based on our previous works|lCJ. Since the SUSY 
phases in the mSUGRA model are so small as to be ignored safely and the Yukawa cou- 
plings are the main source of the flavor structure, we found that the SUSY contributions 
to the CP asymmetries in charmless B decays are generally small except in most B — > PV 
channels. Currently, people have tried to measure the CP violation in B decays, but the 
data are always error- weighted and consistent with zero except for B — > 7r ± fC =F . More- 
over, large theoretical uncertainties also exist for example in the QCDF approach per se 
where the weak annihilations and other potential power corrections are not well calculated. 
Therefore, low experimental statistics and large theoretical uncertainties together prevent 
us from testing the mSUGRA model through studies of the CP-violating asymmetries at 
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present. We are waiting for great progress in both the theory and experiment. 
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